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1. Introduction 



NONLINEAR SCHRODINGER EQUATION ON GRAPHS: 
RECENT RESULTS AND OPEN PROBLEMS 

DIEGO NOJA 

Abstract. In the present paper an introduction to the new subject of nonhnear dispersive hamil- 
tonian equations on graphs is given. The focus is on recently estabhshed properties of solutions in the 
case of nonlinear Schrodinger equation. Special consideration is given to existence and behaviour of 
solitary solutions. Two subjects are discussed in some detail concerning NLS equation on a star graph: 
the standing waves of NLS equation on a graph with a S interaction at the vertex; the scattering of 
fast solitons through an Y-junction in the cubic case. The emphasis is on description of concepts 
and results and on physical context, without reporting detailed proofs; some perspectives and more 
C^ ambitious open problems are discussed. 
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In the last decades a large amount of work has been done concerning existence and behaviour 
^ of solutions of nonlinear dispersive equations of Hamiltonian type. This is in part a consequence 
^ of the fact that many fundamental physical models belong to this family. In particular nonlinear 
^ Klein-Gordon equations and their relatives are milestones of classical and quantum field theory; but 
d in fact an important boost to these developments comes from more phenomenological models, such as 
^ Korteweg-deVries equation describing shallow water waves in certain approximations or the ubiqui- 
tous nonlinear Schrodinger equation which describes electromagnetic pulse propagation in nonlinear 
CN (Kerr) media, Langmuir plasma waves, or in the quantum realm Bose-Einstein condensates, where it 
^ is better known under the name of Gross-Pitaevskii equation. All the above Hamiltonian equations 
QQ share a common characteristic: they admit solitary solutions, or briefly solitons. As it is well known, 
^""^ solitons are solutions emerging due to a balance between nonlinearity and dispersion and they are 
related to symmetries of the equations; in some relevant examples (as one dimensional KdV and 
cubic nonlinear Schrodinger equation) they also have a strong relation to the complete integrability 
C^ of the infinite dimensional Hamiltonian system to which they refer, but their existence and main 
^ interesting properties are by no means restricted to integrable equations. Solitons are nondispersive 
K*- and to some extent particle-like solutions of certain PDEs; similarly to the equilibrium points of 
k> finite dimensional dynamical systems, they are an essential point of departure in the description of 
5-j the phase portrait of the model equation which they solve. In the present paper it will be shown the 
existence and role of solitons for some model PDEs on ramified structures, in particular the simplest 
kind of them, the so called star graphs. The model equation described in some detail is the nonlinear 
Schrodinger equation on graphs, about which a rigorous mathematical activity is developing. Before 
describing mathematical models let us give a look at the literature of physical origin which motivates 
the study. Two main fields where NLS equation enters as a preferred model are optics of nonlinear 
Kerr media and dynamics of Bose Einstein condensates. Both of these quite different physical situa- 
tions have potential or actual application to graphlike structures. In nonlinear optics we can mention 
arrays of planar self focusing waveguides, and propagation in variously shaped fiber optics devices, 
such as Y-junctions, H-junctions and more complex examples can be considered. Papers relevant to 
these items are [U [2] , where also symmetry breaking phenomena related to geometry are studied. In 
[21 H] experimental evidence of interaction of solitons with inhomogeneities and defects is given, in 



particular scattering and capture of solitons in photonic traps and escaping of solitons from poten- 
tial wells. These last phenomena are studied on one dimensional media but they suggest a natural 
generalization to simple graphlike structures, such as Y-junctions. In [5], after an analysis of general 
issues discussed later in this paper, an example of potential application to signal amplification in 
resonant scattering on networks of optical fibers is given. In the field of Bose-Einstein condensates 
and more generally of nonlinear guidance of matter waves there has been an increasing interest in 
one dimensional or graphlike structure. Boson liquids or condensates can be treated in the presence 
of junctions and defects in analogy with the Tomonaga-Luttinger fermionic liquid theory, with appli- 
cations to boson Andreev-like reflection, beam splitter or ring interferometers (see for example [HI [7j 
and references therein). Apart experimental activity, some other theoretical and numerical studies 
should be mentioned. NLS equation on graphs has discrete analogues (see [H [9l [TOl E]), where 
several dynamical behaviours can be studied analytically and numerically. Other discrete models are 
spin models on star graphs related to Kondo model (see [I2])- Possible integrability of cubic NLS on 
star graphs with special boundary conditions is discussed in [T3]) and finally quantum field theory on 
star graphs with a consideration of symmetry, integrability and analysis of special models is studied 
in [H]. 

1.1. Preliminaries and the mathematical model. To pose the models, we briefly recall some 
preliminary notions (see[T5l [161 HZl UHl UHl l2Qj for a more systematic account). By a metric graph Q 
it is meant a set of edges {ej}^^^ and vertices, with a metric structure on any edge. Every edge of 
the graph is identifyied with a (bounded or unbounded) oriented segment, Cj ^ Ij- A function on a 
graph is a vector 

\I' = (ipi, ...,iPn) with ipj = ipj{xj) ; Xj G /-,■ . 
In the following, as above, we will denote the elements of L'^{Q) by capital Greek letters, while 
functions in L^(]R+) are denoted by lowercase Greek letters. The spaces L'^{Q) are defined in the 
natural way by 

N / N 

Only for the L^-norm we drop the subscript and simply write || • ||. Accordingly, we denote by (■, ■) 
the scalar product in L'^{Q). When an element of L'^{Q) evolves in time, we use in notation the 
subscript t: for instance, "^t- Sometimes we shall write \E'(t) in order to emphasize the dependence 
on time, or whenever such a notation is more understandable. In a similar way one can define the 
Sobolev spaces 

N 

i=i 
The definition however should be used with caution, because as such, without specification of be- 
haviour at the vertices, these spaces have not the usual properties of Sobolev spaces, in particular 
continuous or compact embeddings. The main example of metric graph that will be considered 
here is a star graph, which is characterized by a single vertex v and N infinite edges, and we put 
Cj ~ (0, +cxo), V = 0. 

After functional spaces, differential operators can be given on graphs. In particular we are interested 
in operators connected with the Laplacian and some variants. On the Hilbert Space 

N 

M = 0L2((O,oo)) 
i=i 
with elements 

^ = (V'i,...,^7v)eH 
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endowed by a norm 
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we consider the operator Hg 

Ha^ - . 
on a suitable domain 

^(ifg) = ^RiJ^((0, oo)) & self-adjoint conditions at the vertex 
i=i 
The choice of the boundary condition of course qualifies the operator, and different boundary con- 
dition give rise to different dynamics when the operator Hg is the generator of an evolution, e.g. 
Schrodinger or Klein-Gordon. 

As it is well known, unitary N x N matrices U parametrize the family of selfadjoint Laplacians 
on Q, and the relation between U and the specific boundary condition at the vertex is given by 

/^i(0)\ /^'i(0)\ 

(f/-l) : +^(f/ + l) : =0 

However we are not interested here in the greatest generality, and we restrict ourselves to the 
following simple cases: 

notice that for a graph with two edges, i.e. the line, continuity of wavefunction and its derivative 
for an element of the domain makes the interaction disappear; this fact justifies the name of free 
Hamiltonian. 

- S condition: [Ujk = 2{N + ia)~^ — 6jk] 

N 

^{v)=^i{0) = MO) = --- = ^n{0), ^7/';.(0)=a7/'(0) aeR 

i=i 
The S condition includes the so called Kirchhoff condition a = as a special case, but it is convenient 
to distinguish between them. Notice that for a graph with two edges, i.e. the line, continuity of 
wavefunction and its derivative for an element of the domain makes the interaction disappear; this 
fact justifies the name of free Hamiltonian. 

From now on Ha is the operator Hg with delta condition in the vertex with "strength" a. To simplify 
notation the Kirchhoff case will be indicated with the symbol H. Another more singular interaction 
to which we will refer occasionally is given by the 

- 6g condition: 

5^^;.(0) = 0, z^,(0)-^,.(0) = ^(V^;.(0)-^UO)), j,/c = l,2,...,n, 
i=i 

which coincides, in the case of the line, with a S' interaction of strength /3/2. The operator H^ 
has simple spectral characteristics. The absolutely continuous spectrum coincides with [0, oo) . As 
regards the point spectrum, 

- For a < 0, CpiH^) = -a^N^ 

- For a > 0, aplnj) = {0} 

The presence/absence of the (negative) eigenvalue is the motivation for the name of attractive/repulsive 
delta vertex given to Ha in correspondence to the cases a < / a > 0. More specifically, a 6 vertex 



with a < can be interpreted as modeling an attractive potential well or attractive impurity. In 
fact, as in the case of the line, the operator H^ is a norm resolvent limit for e vanishing of a scaled 
Hamiltonian H^ = H + aV^, where K = -V{-) and l^ is a positive normalized potential defined on 
the graph in the natural way and H is the free Hamiltonian ([17j and reference therein). 
Finally, the quadratic form associated to Ha is 

Q[^] = ^ii^'ir+f 1^(0)1 

^{Q) = {^eH' s.t. ^(0) = ^i(O) = ... = tjjNm = S 
and 

Notice that 3l[Q\ the form domain of Q, is independent on a. In the following we will use the 
convenient notation £^, calling it the energy domain (it is often considered the "true" Sobolev space 
of order 1 on ^). 

To introduce the nonlinearity, we define a vector field G = (Gi, ■ ■ ■ ,Gn) '■ C" — t- C" acting 
"componentwise" as 

^.(C) = ^(101)0 with (7 : M+ ^ M and C = (C.) e C" . 

The vector field G enjoys the important property of gauge (U(l)) invariance, i.e. G{e'^^Q = e^^G{Q . 
After this preparation the more common evolution equations on the graph can be defined in the 
obvious way. Main examples are 
i) the nonlinear Schrodinger equation 

(1.1) i — ^t = i^a^i + G'(^t) (NLS equation) ; 
ii) the nonlinear Klein-Gordon quation 

(1.2) - —^t = Ha^t + m'^'t + G{^t) (NLKG equation) . 

Notice that from a mathematical point of view the nature of a PDE on a graph amounts to a 
system of PDE's on suitable finite or infinite intervals (in the case of a star graph N halflines) in 
which the coupling is given exclusively through the boundary conditions at the vertices. In the 
following paragraph we will see several illustrations of this simple remark in the case of our main 
example, the NLS equation. 

A further specification of the model is given specializing the study to the important case of a power 
nonlinearity 

g{z) = ±|^|2^ /i > . 

In the case of NLS equation, the minus sign corresponds to so called focusing nonlinearity and the 
plus sign to the defocusing one, both meaningful in the applications. In the NLKG equation, the 
minus sign is the most relevant. The cubic case is especially important both for the NLS, where 
it describes the most common phenomenological situations in nonlinear optics and BEG, and in 
the case of NLKG equation, where it corresponds to the quartic interaction in field theory. Simple 
combination of monomial nonlinearities, such as the physically relevant " cubic-quintic" nonlinearity 

g{z) = — 1-2|^ + e|2;|^ 

encountered in nonlinear optics and other fields (see for example [21]) can be considered as well. 
In the following we almost exclusively refer to the example of NLS with focusing power nonlinearity 
and an attractive 6 vertex. After fixing the model the first preliminar and essential information 
regards its well posedness, i.e. existence (local or global) and uniqueness of solution in suitable 



functional spaces. For our preferential model, the NLS equation, the classical line of attack to well 
posedness is through the integral form of the equation, given by 

ft 



;i.3) ^t = e-'^"*^o + ^ f e-'^-^'-'^l'^sl^^'^sds 

Jo 



A formal analysis show that for NLS on graphs there exist conserved quantities, similarly to the case 
of the line or of open set in M". These are the mass 

(1.4) M[^] = ll^f ; 

and the energy 

(1-5) Em = l\mr - ^11^112^2 + f k(o)r ; 

this last quantity for fi G (0, 2) is well defined on the domain ^{E) = £ oi the quadratic part of the 
energy, i.e. the energy of the linear Schrodinger dynamics on the graph. 

A sample rigorous result giving well posedness of NLS equation on star graphs and conservation 
of mass and energy is the following. 

Theorem 1 (Local and global well-posedness in £). For any \l'o G £, there exists T > such that 



the equation (1.3) has a unique solution ^ G C''([0,T),£) n C^([0,T),£^* 



Moreover, eq. (1.3) has a maximal solution xj/™**^ defined on an interval of the form [0,T*), and the 



following "blow-up alternative" holds: either T* = oo or 



lim ll^r^'IU = +00, 



where we denoted by vl/^a-^ the function v]/™^-^^ evaluated at time t. 

Moreover, in the same hypotheses, the following conservation laws hold at any time t: 

M[^i] =M[^o], E[^t] = E[^o] ■ 

Finally, for < fi < 2 and any \E'o G £ , the equation ( |1.3[ ) has a unique solution \E' G C°([0, oo), £^) fl 
C'{[0,oo),£*). 

In the previous theorem, £* is the dual space of £" . As in the case of NLS equation on the line, 
we will call the range of nonlinearities < /i < 2, where existence for all times is guaranteed, 
subcritical case, in contrast to the supercritical nonlinearities (/i > 2) where blow-up occurs (in the 
form described in the theorem), or the threshold critical (/i = 2) nonlinearity where global existence 
depends on the size of the initial datum. For details on the proof and generalization to more general 
couplings at the vertex of the star graph see [221 123]. The same theorem holds true for more general 
nonlinearities with essentially the same proof. 

The previous examples are modeled starting from generators which are nonlinear perturbations of 
" quantum graphs" . 

On the other hand, in principle, other equations could be considered. Without embarking here in 
a general theory we limit to mention some special models. The first one is given by the Benjamin- 
Bona-Mahony equation (BBM) describing the unidirectional shallow water flow under the long wave 
and small amplitude approximation, 

du du du d^u 

dt dx dx dx'^dt 
Of course here the unknown u is real. It is immediate to formally extend such an equation on a 
graph. To simplify the exposition, consider the case of a star graph with three edges, an Y-junction. 
Let us define a vector u = (mi,M2,M3) where Ui : (0, oo) — )■ M and suppose that at the vertex the 
three components Ui satisfy Kirchhoff boundary conditions. These conditions are rather natural 



in the context of water (and other fluids) waves, corresponding to continuity of the flow and flux 

balance. 

Now orienting the edges of the Y-junction in the outgoing direction as above and setting cti = 

— 1 , (72 = (Js = 1, let us consider the evolution problem 

dui dui dui d^Ui ^ , 

ot ox ox ox^ot 

u, e H\R^) , uM = um = u-M , ^(0) + ^(0) + ^(0) = . 

ox ox ox 

This is a system of scalar BBM equations on the halfline, coupled through the Kirchhoff boundary 
condition at the origin. In [2^ a generalization of this systems to trees (including star graphs) is 
studied as a model of cardiovascular system, and in particular it is shown that BBM on a tree is 
well posed. It is interesting that (see |2S]) travelling waves for BBM equation on a tree have been 
recently constructed for particular vertex conditions. By the way a result of well posedness of NLS 
on trees in not yet proven, although well plausible; a step in this direction with proof of relevant 
dispersive estimates is given in [261 EZ]- Other nonlinear models of quite different nature are the 
reaction-diffusion equations on networks, about which some literature exists, in particular regarding 
pulse propagation in axons and neural networks according to the FitzHugh-Nagumo model and its 
variants (see [28] and references therein). Concerning the mathematical setting, a last word should 
be said about the one dimensional approximation given by the NLS equation on a graph. The graph 
should be a limit in some sense of a more realistic systems of thin tubes (or guides) connecting at 
junctions. A first problem is getting the limit of a certain dynamical model, e.g. NLS equation defined 
on the system of thin guides when transversal size of tube vanishes. It is reasonable to conjecture 
that this limit should be a NLS on a graph as defined before, but with which boundary conditions? 
And how the boundary conditions at the vertex could depend on the limit process of shrinking the 
tubes, guides and junction size? A further difficulty could be the dependence of limit process on the 
conditions at the boundary of the thin tube, e.g. Neumann or Dirichlet. These problems have been 
tackled, with partial solution, in the linear case (see [291 [301 [31] ^ind reference therein). They remain 
open for nonlinear models, where there is no literature, with the only exception (to the knowledge of 
the present author) of a series of papers about the reduction of the Ginzburg-Landau equation and 
its stationary counterpart from thin tubes to graphs (see [321 I3SI [31] and references therein), where 
some special boundary conditions at the junction appears in the limit. A second point of view is less 
phenomenological, and related to the deduction of evolution equation for Bose Einstein condensates 
from first principle. The dimensional reduction of BEC using scaling trapping potentials is a well 
understood process in several limiting regimes; in particular a Gross-Pitaevskii energy functional on 
the line describes the so called cigar-like BECs under certain conditions (see [351 [M| and reference 
therein). Being confident that a Gross-Pitaevskii equation realizes the correct quasi one dimensional 
limit, a similar procedure could be attempted on graphlike structure, for example on a Y-junction; 
to simplify the analysis, one could start directly from a A^-body theory on the graph, and to attempt 
at taking the A^ — )■ oo limit with suitable scalings. Notice that the main problem, i.e. the treatment 
of the boundary condition, it is open also in the simplest case of a graph with two edges, i.e. a line 
with a defect. 

2. The NLS equation on star graphs: rigorous results 

After setting the mathematical model and giving the main physical premises and possible applica- 
tions, we turn to a description of some of its dynamical features. Again we refer to the case of the 
focusing NLS equation where more information it is known. Two main topics have been rigorously 
indagated in the last years: the existence and characterization of standing waves, and the scatter- 
ing of solitons through a juctions. To introduce the two subjects, let us preliminarly recall some 



properties of NLS solitons on the line. Let us consider the usual focusing NLS equation with power 
nonlinearity in one space dimension 

dzt 11 

(2.1) i — {x,t) = -—{x,t)-\u{x,t)\^f'u{x,t) xeM,t>0; 

as it is well known, it admits a special solution u{x,t) = e**0 with 

(2.2) (f){x) = [(/i + l)]57i sechM(/ix) . 

A richer family of solution is obtained through application of Galileian and scaling symmetries of 
the NLS equation: 

(2.3) u,,,,,^{x, t) := e'5"e-**'^e'"*a;5^0(v^(x - xq - vt)) . 

Notice that frequency uj of the oscillation and amplitude of the solitary wave Ux^^v^u) are nonlinearly 
related, and in particular the greater is the amplitude, the greater is the oscillation frequency. 

2.1. Standing waves of NLS on star graphs. Let us begin with the first subject, standing waves. 
On the line they appear putting f = in the previous family of solitary solutions, and they have 
the character of localized solutions (or "pinned" solitons, in the physical literature) around a certain 
centre Xq. In particular they are the only soliton solutions when travelling waves are excluded by the 
presence of inhomogeneities; these can be represented for example by external potentials, magnetic 
fields or, as in our case, by a boundary condition at the junction. In these cases the localization of 
the standing wave is around stationary points of external potentials, or at the location of singular 
interactions in case of junctions, point defects, etc. Here we will define standing waves as finite 
energy solution to a NLS equation (for other models with f/(l) symmetry, such as NLKG equation, 
the definition is exactly the same) of the form 

(2.4) mt{x) = e*"* ^^{x) . 

The function $^ is the amplitude of the standing wave. In particular we are interested in standing 



waves of the NLS equation on graph, equation (1.1). A regularity argument shows that the standing 



waves belong in fact to the operator domain of Ha- Correspondingly there should exist a frequency 
u and an amplitude $^ which satisfy in the strong sense the stationary equation 

(2.5) i7,$<. - |<I>J2^<I>^ = -0;$^ . 

The analysis of this equation on a star graph is simple (see [23l |37] for details). 
On every edge the operator H^ coincides with the second derivative, and so it holds 

the most general solution with cf) G L^(]R"*") is (introducing explicitely the dependence on parameters) 

1 1 

(2.6) 0(o", a]x) = a [(/i + l)u\ '^>^ sech'' {^y/uj{x — a)) , |(t| = 1 a G M 

so that the components of the amplitudes are 

(*c^)i = (t){(Ji,ai) , 

where (Jj,aj have to be chosen to satisfy $(^ G V{Ha) ■ Notice that on every edge i the stationary 
state has an amplitude which is "bump-like" or "tail-like" in shape, according to the position of the 
center Oj , within the edge i or not. 
Continuity at the vertex implies the following conditions on the parameters of the amplitude 

(Tj = 1 , aj = EjQ , Ej = ±1 j = 1,- ■ ■ N , a >0 . 

Being aj = 1 Vj we will drop its mention in the following. 







Figure 1. A^ = 3 stationary states. On the left a < 0; on the right a > 0. 



Moreover, imposing the 6 vertex boundary condition to $0; gives 



N 



(2.7) 



tanh(/iA/cJa) >^ Ei 



i=l 



a 



so that ^1,1=1 ^i must have the same sign of a if a 7^ 0. 
Immediate quahtative consequences of the above hmitations are: 

- a > strictly more bumps than tails 

- a < strictly more tails than bumps 

- a = same number of tails and bumps or a = 

- For any value of a 7^ there are [^^^] states ([s] is the integer part of s). 

2 

- Lower bound on the allowed frequencies: ^ < u 

We index the stationary states ${, with the number j of bumps. With the above limitations this 
identifies completely the state. More explicitely the j-bumps state <J>^ is given by 



(2.8) 



{K)S^) 



(a-^; x) 
(—a^\ x) 



1,...J 

j + l,...iV 



(2.9) 



1 



arctanh 



/iA/W 



a 



{2j-N)V^ 



Concluding, solutions of ( |2.5[ ) for a > are given by $^ with j = [N/2 + 1] , . . . , A^ and for a < by 
$iwithj = 0,...,[(iV-l)/2]. 

The situation is pictured in Figure 1 for the A^ = 3 star graph. 

So for a 7^ , for every N and w > -^ there exist branches {^i,} of stationary states (the branch 
is unique only in the case N = 2, i.e. the line). More precisely, the state with j = arises for 
while to have states with j > higher frequencies are needed, according to the general 



u > 



Ar2, 



relation u > jj^^^- 

Now let us consider the Kirchhoff case, a 



0. 



From an analysis of the boundary conditions it follows that star graphs with odd or even number of 
edges behave differently. 





Figure 2. The Kirchhoff case a = for odd (left) and even (right) number of edges. 

For N odd the only value of a compatible with the boundary conditions is a = 0. So the stationary 
state is unique 

($°)^(x)=0(O,x) z = l,...,N 

and it it is composed by N half solitons continuously joined at the vertex. 

For A^ even every real value of a is compatible with the boundary conditions and there is the same 
number of tails and bumps on the graph. A one-parameter family of stationary states exists and is 
given by (notice the slight change of notation used for the Kirchhoff solitary wave only) 



i^: 



[X 



[-a,x) i = l,...N/2 
[+a,x) i = N/2 + 1,... N 



a G 



These stationary states could be thought as N/2 identical solitons on N/2 lines. The situation is 
depicted in Figure 2. As a consequence there exist travelling waves on a Kirchhoff graph with an 
even number of edges: simply translate the complete solitons on every fictitious line by the same 
amount vt: 



Notice that in the Kirchhoff case stationary states exist for every positive u. 

A similar construction can be performed for NLKG standing waves with different restrictions on 

parameters, i.e. \uj\ < \ m? — -^ . Details will be given elsewhere. 



2.2. Variational properties of standing waves. After constructing the stationary states, the 
natural problem is to identify the ground state and, if possible, to order the states in energy, i.e. to 
describe the nonlinear spectrum of the NLS equation on the graph. This is a variational problem, 
which is also relevant for the analysis of stability of standing waves. A difficulty immediately arises. 



in that the NLS energy (1.5) for the focusing NLS equation is unbounded from below, as easily 



recognized (this is not the case for the defocusing nonlinearity) . So the seemingly natural definition 
of the ground state as the minimizer of the energy is meaningless. Nevertheless the physics of 
the problems behind the NLS equation suggests that a possible relevant variational problem is a 
constrained variational problem: to minimize the energy at fixed mass. With this constraint and for 
subcritical nonlinearities fi < 2 the energy is shown bounded from below for every finite energy state, 
as it is in the case of NLS equation on M". In fact the following result holds true (see [38] for details 
and proofs) for the focusing NLS on a star graph with an attractive delta vertex, a < 0. 

Theorem 2 (Minimizers for the Energy functional). Let m* be defined by 



(2.10) 



m 






1/a* 



\a\ 



(1-t^ 



--1 



dt . 
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Let a < and assume m < m* z/0</i<2 and m < min{m*, '^^^^ } if n = 2 and set 

-V = inf{E[$] s.t. $ G ^ , M[$] = m} . 

Then < z/ < oo and there exists $ such that M[<^] = m and E[^] = —v. 

Moreover the minimizer $ coincides with the N-tail state $°^ where Uq is such that M[$[J,jJ = m. 

So, for every mass above a certain threshold m* (which however it is not optimal) the problem 
of minimizing the NLS energy on the graph at constant mass has a solution if the vertex carries 
an attractive S interaction. More precisely, if the mass constraint coincides with the mass of A^-tail 
state, the minimizer is exactly the A^-tail state. Some comments are in order. 
After some calculation the mass of the stationary states as a function of u turns out to be 



(2.11) M[$i] = il^±^J-^ 



(iV-2j) [ {l-t^)-^-'dt + 2jI 



(JV-2j)c^2 



where / = /(/i) is a certain constant depending on fi only. We recall that ${, is defined for u G 



2 



{N-2jr ' ' 



oo j. So from (2.11) one easily concludes that the functions M[$^] are increasing in u and 



lap 



the minimum value is in correspondence to the threshold (j^_2-)2 - As a consequence the iV-tail state 
can have an arbitrarily small mass while the other stationary states have a minimal mass separated 
from zero. Stated otherways, the manifold M[\E'] = m for m < m* may not contain all the stationary 
states, due to the fact that some of them could have too large masses. On the contrary the A^-tail 
state always belongs to the constraint manifold since its mass has vanishing lower bound. Taking 
into account the dependence of m* on a, one concludes that for small |a| the constraint manifold 
contains only the N-tail state while for large |a| all the stationary states belong to the constraint 
manifold. Moreover from the expression of m* it follows that to guarantee the existence of ground 
state for a given mass constraint one has to have a sufficiently deep 6 well. So alternative statements 
and proofs of the above theorems are obtained fixing m and requiring a to be sufficiently negative. 
Analogous remarks also apply to the critical case /U = 2, the quintic NLS. 

With the above precisations a well defined order in energy exists for the stationary states: the energies 
of the stationary states $fj, where u = Uj is such that M[<l>^,] = m, are increasing in j, i.e. they can 
be ordered in the number of the bumps (see [38] . 
At least in one case things are simple: the cubic case. If /i = 1, then Uj is independent of j and 

, (m + 2|«|)2 

^^■"^ =Sa^' 

So, in the cubic case the energy spectrum at fixed mass can be explicitly computed: 

N s 1 lap 1 (m + 2\a\f 1 laP 



3 3 (2j - iV)2 2A m 3 (2j - N)^ ' 

The energy of the ground state is given by 

ek*] = -^^(^' + 6^i«i + i2i«r) . 

As a final remark, notice that the ground state of the system is the only stationary state which is 
symmetrical with respect to permutation of edges. 

A second variational problem which has both a physical and mathematical relevance is related to the 
minimization of the action functional. The action for a nonlinear Schrodinger equation is obtained 
adding the non linear term to the usual action of the linear Schrodinger equation 

'•t2 



A[^] = r (^ f gm(^ ^) dx + E{^U dt 
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This expression, for a standing wave "^t = e^^*(^^ takes the form 

^[e^"*$J = (ti-to)(i5;[$.] + |||$J|^) . 
This fact suggests to consider the "reduced" action 

(2.12) s^m = E[^ + |m[$] = \\wr + |ii«i>f - ^;^ii*ii2;::2 + f i^wp , 

which we continue to call action with a common abuse. 

Apart from this lagrangian origin the action S^^ has, at least in the contest of BEC, the physical 

interpretation of the the grand-potential functional of the condensate corresponding to the chemical 

potential u. 

Whatever the theoretical interpretation, the action functional just introduced enjoys the important 



property that its stationarity condition 5'^[$] = coincides with the stationary equation (2.5). As 



for the energy it is easy to see that the action is unbounded from below. Nevertheless a solution of 



equation (2.5) satisfies necessarily (just take the scalar product of the stationarity equation with $) 



to the constraint 



m = \w\? - wntxi + ^m? + (^\m? = s'^^ = o 



\2^l+2 

So it is expedient to search for minima of the action restricted to the above natural constraint, 
also called in mathematical literature Nehari manifold. It contains all the stationary states by very 
definition. Now an immediate calculation show that restricted on the natural constraint the action 
is 

Q [$1 - _^^||$||2m+2 

^"^*J"2/i + 2"*"2^+2' 
and so it is bounded from below and nonnegative (this is true for every power nonlinearity /i, notice 
the difference with the constrained energy). The absolute minimum of the action constrained to the 
Nehari manifold, if existing, is called ground state of the action, as for the energy. Notice that by 
Lagrange multiplier theorem, the stationary points of the energy at fixed mass are stationary points 
of the action with Lagrange multiplier uj. In fact the ground states of the two problems coincide, as 
a consequence of the following result (see [23] for a complete discussion and proof). 

Theorem 3 (Minimizers for the Action functional). Let ^ > 0. There exists a* < such that for 
—Ny/uJ < a < a* the action functional S^^ constrained to the Nehari manifold admits an absolute 
minimum, i.e. a $ 7^ such that /a;[$] = and S^^l^] = ini {Suj[^] : Ilu[^] = 0}. Moreover, for 
—Ny/u < a < a* the ground state is <I> = $° . 

The threshold a* in the above theorem is known as a function of N, ^ and a. The proof of Theorem 
|3] is quite different from that of Theorem [2j Nevertheless the origin of thresholds a* and m* is 
analogous, and it relies on the fact that the action of NLS on a Kirchhoff junction has no ground 
state: the infimum exists, but is not attained at any finite energy state. To explain, let us consider a 
Y junction, i.e. a A^ = 3 star graph. There exist sequences of "runaway" states (see Figure 3) on the 
Kirchhoff graph given by a complete soliton on a couple of halflines plus a correctly joined soliton 
tail on the third halfiline; Kirchhoff boundary condition are easily verified, so we have a sequence of 
domain elements. Now, the more the big soliton shifts to infinity and the tail estinguishes, the less 
is the action, which can be shown to converge toward its infimum, strictly lower than the action of 
any stationary state. The same phenomenon occurs at the constrained energy level as shown in [3U] . 
The bad behaviour of the Kirchhoff junction {a = 0) prevents the action with a small, or the energy 
with a big m, to have a constrained absolute minimum. It is a conjecture that the action has a local 
constrained minimum that is larger than the infimum when the condition —Ny/u < a < a* fails, 
and similarly for the constrained energy with a corresponding condition on mass. 
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Figure 3. A runaway state on the graph 

2.3. Orbital stability of standing waves. Stability is an important requisite of a standing wave, 
because unstable states are rapidly dominated by dispersion, drift or blow-up, depending on dynam- 
ics, and so are undetectable (stability and instability of NLS with a 6 potential on the line is studied, 
partly numerically, in |10]; see also reference therein). Due to gauge or U{1) invariance of the action 
and dynamics, a standing wave is not stable in the usual Lyapounov sense. This is a general fact 
in the presence of symmetries, and it is well known in the example of relative equilibria for finite 
dimensional mechanical systems. In our case to introduce the main concept af orbital stability, let 



us consider the special solution e 0(x) to equation 2.1, where is the initial datum given in 2.2 

by scaling and U{1) invariance, u{x,t) = e'^^*uJ~i^(|){^JuJx) is the solution corresponding to the initial 

datum ujT^ (f){^/uJx) . Choosing uj close to 1 the two initial data are close. But their time evolutions 
are not, because of different frequencies which make distance of solutions to vary with time: 

sup I |u(x, t) - e'^(t){x) \\m > 1 101 1^1 • 

The same phenomenon occurs on M for Galileian invariance: slightly different velocities make trav- 
elling solitary waves separate each other. In the case of graphs Galileian invariance is not relevant 
because it is broken by the junction, travelling waves do not exist and so we concentrate on standing 
waves. The point of the previous discussion is that the stability has to be defined up to symmetries: 
solutions remain close to the orbit 6*^$^^ of a ground state for all times if they start close enough to 
it (see for a general discussion [HI US]). The orbit of $a; is 

Definition 2.1. The state $(^ is orhitally stable if for every e > there exists 5 > such that 

{!(«*„))<« ^ W).^(*„))<. w>o 

where 

d(^,^($,))= inf ||^-m||^. 

The stationary state <l>a; is orbitally unstable if it is not stable. 

A general theory of orbital stability was established in the eighties in the classical papers by 
M.Weinstein j43j and Grillakis-Shatah-Strauss [^ 1^ . and developed in a number of subsequent 
papers by many authors; it applies to infinite dimensional Hamiltonian systems (such as abstract 
NLS equation) when a regular branch of standing waves uj i— )■ \E'^, not necessarily ground states, 
exists (see for example jlHl ST] for recent surveys and results about stability and instability in NLS 
case). The first step is to give the NLS equation on graph an Hamiltonian structure. This is achieved 
in a standard way, considering an element of L^(^, C) as the couple of its real and imaginary part. 
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\i/ = U + iV = {U,V) and endowing the Hilbert space L^{g, C) = L'^{g, M) © L'^{g, R) so obtained 
with the real scalar product (\E'i, "^2) = 5^^ L \E'i\E'2 = fg U1U2 + V1V2, and analogously decomposing 
the higher Sobolev spaces. The NLS on a graph turns out to be a Hamiltonian system 

(2.13) i,{v)-^my\' ^=(-Vo)- 

where E[U, V] = E[^] and the functional derivative is defined as usual 

E'[U, V]{X, Y) = jE[{U, V) + e(X, r)],=o V(X, Y) e H\g, R) © H\g, R) . 

Now linearize the Hamiltonian system around the ground state setting 

(2.14) ^{t)= {^l + W + iZ)e''^' . 

Notice that the previous definition amounts to pass to a rotating coordinate system, comoving with 
the ground state. Then the fiuctuations W and Z satisfy 

where 



^' Z J \C2Z 

and £1 and £2 are matrix s.a. operators: their domain coincide, with a slight abuse of notation, 
with ^{Ha) and the action is given by 

(2.16) (A),, = (^-^ + c - (2/i + 1)|<,P'^) 5,, , 



Due to the definition 2.14 , it turns out that the linearized operator C coincides with the second 
derivative of the action: C = S'^($°), after identification of the sesquilinear form with the operator 
via the scalar product. 

The first information about stability is given by the spectrum of the operator C. When the linearized 
operator C admits at least one eigenvalue with nonvanishing real part the stationary state $J^ is said 
spectrally unstable, otherwise it is said spectrally stable. In fact, due to the presence of conservation 
laws, in particular of the mass H^E'lP, the system can be spectrally unstable without being orbitally 
unstable. Precisely, according to the general theory of Weinstein and Grillakis-Shatah-Strauss, for a 
solitary solution <l>^ being orbitally stable is sufficient that: 

i) spectral conditions hold: 

ii) ker £2 = {^uj} a^nd the remaining part of the spectrum is positive 

22) ker£i = {0} 

23) the number of negative eigenvalues (the Morse index) of £1 is equal to 1 
ii) Vakhitov-Kolokolov condition -^jW^uiW^ > holds. 

These conditions can be verified and one obtains the following theorem 

Theorem 4. Let 0<fi<2,a<a*<0,u>^. Then the ground state $° is orbitally stable in 

E . Moreover, if ^ > 2 there exists u* such that $° is orbitally stable for u G (^,aj*) and orbitally 
unstable for u > u* . 
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The theorem gives orbital stabihty of the ground state for every u also for the critical nonlinearity 
fi = 2 . The case u = u* is undecided. The proof is a calculation for the Vakhitov-Kolokolov 
condition; and concerning the spectral conditions, the Morse index of £i is one as a consequence of 
the fact that the action has a minimum at <l>° restricted to the codimension one Nehari manifold. 
The other details are in [23] . 

We end this subsection with some comments. There is a second strategy to show orbital stability, 
which again makes use of a variational property. A stationary state which minimizes the energy at 
constant mass is orbitally stable (see the classical paper |12], where several example are treated). So a 
direct consequence of Theorem [2] is orbital stability of ground states. Nevertheless some remarks are 
in order. The first is that while Theorem [2] and in general the concentration compactness technique 
developed in |12] give information only if absolute constrained energy minima (i.e. ground states) 
exist, the above Weinstein and Grillakis-Shatah-Strauss theory is more general, and it is in principle 
applicable to every stationary state of the action, for example excited states of NLS equation on star 
graphs above described. For excited states the expectation is orbital instability, which is in fact the 
case, as it will be shown elsewhere. The difficult part of the analsis is the calculation of the Morse 
index of operator Ci and the use of general results in [IS] and their recent refinements in [17] . One 
could ask if a similar analysis could be performed on more complex graphs, for example trees or 
graphs with loops. Of course the fact that stationary states on star graphs are completely known, 
which is a rare case, is a strong facilitation in obtaining precise results. In the case of less trivial 
graphs it is generally impossible to obtain explicitly standing waves, but some simple non trivial 
graphs can be probably treated along the lines discussed. On more general grounds, when the linear 
part of the model, i.e. the underlying quantum graph has an eigenvalue, for example corresponding to 
the linear ground state, bifurcation theory suggests that a branch of nonlinear standing waves exists 
and it bifurcates from the vanishing solution in the direction of the linear ground state (see [38] for 
a classical application to the case of NLS with external potential). But there are several problems 
which arise at this point. The first is that a direct analysis of the conditions guaranteeing orbital 
stabihty or instability of standing waves, in particular counting of negative eigenvalues of £i and 
verification of Vakhitov-Kolokolov condition, becomes in general impossible or at least very difficult. 
A second problem is that bifurcation theory allows to identify branches of nonlinear stationary states 
which have a linear counterpart, but how to obtain branches of states without linear counterpart, 
which exist as we know from the example of nonlinear excited states in star graphs? A guess 
is that excited or in general bound states without linear counterpart bifurcate from (not small) 
solitary waves turning on the external potential. A different extension could be in the direction of 
different boundary conditions at the vertices. In such a case one expects new dynamical effects, for 
example bifurcation and symmetry breaking of ground states as proved for the NLS on the line with 
6' interaction (see [32]). A final open and difficult problem is the so called asymptotic stability of 
standing waves. In our case a standing wave "I'^e*'^* is said asymptotically stable if for every solution 
u{t) starting near ^^ in the energy norm, one has the representation 

(2.18) uit) = e*"-*$^^ +Ut*^^ + R^t) , 

where Ut is the unitary evolution of the linear s.a. operator H^, and \E'oo, Roo{t) ^ L'^iQ), with 
||-Roo(^)|| = 0(t~^) as t — )■ +00 , for some /3 > and u^o > fp ■ So every solution starting near an 
asymptotically stable standing wave is asymptotically a standing wave (not necessarily the original 
one) up to a remainder which is a sum of a dispersive term (a solution of the linear Schrodinger 
equation) and a tail small in time. The physical interpretation of the concept is that dispersion, or 
radiation at infinity, provide the mechanism of stabilization, or relaxation, toward the asymptotic 
standing wave or more generally solitons. See for example [50] which treats NLS equation with a 
potential on the line, and references therein. The problem is very hard and there is only partial 
information also in the case of a NLS equation on the line with a 6 potential (see [51] ) . Part of the 
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difficulties are due to the fact that the Hamiltonian structure plays a role in the analysis and this 
makes unavoidable to study the spectrum of the not selfadjoint and not skew adjoint operator J C 
(the Hamiltonian linearization) and to get dispersive estimates about its evolution exp(tj7£); this 
introduces some interesting and new mathematical problems about operators on graphs, already at 
the linear level. Moreover the proof of possible asymptotic stability requires a control of the decay of 
nonlinear remainders; this control depends in a critical way on several analytical tools (in particular 
dispersive Strichartz estimates) which at present work only for restricted classes of nonlinearities; for 
example, for subcritical power nonlinearities the procedure fails. 

2.4. Scattering of fast solitons on junctions. In the preceeding section the behaviour of localized 
standing soliton solutions of the NLS equation and of the solutions in their vicinity has been studied. 
Here, concentrating on the case of cubic NLS equation on a three edge star graph, we explorate a 
different region of state space of this model, that of the asymptotically travelling waves. 

As recalled before among the family of solitary solutions of NLS equation on the line given by the 
action of the Galilei group on the elementary function 

(2.19) 0(x) = v/2cosh"^a;, xgM, 
there are the translating waves 

(2.20) 0xo,i,(a;, t) = e*t^e-**Te*V(x - Xq - vt) xeR,teR,veR. 

This special state, when put on a single edge and pushed to infinity, could be reasonably considered 
as an asymptotic soliton traveling on the graph. Of course the presence of the vertex breaks galileian 
invariance and the soliton cannot rigidly translate in the course of evolution, also for the simplest 
graph, i.e. the star graph with Kirchhoff boundary condition at the vertex. The classical and well 
known algebraic and analytic techniques to construct exact solutions of (integrable) cubic NLS on 
the line fail on a graph. The interaction with the Y-junction could be in general quite complex, and 
from a mathematical point of view essentially nothing is known, if not in the special case where the 
asymptotic solitary wave is a fast soliton, in a sense that will be made precise later. In such a case, 
after the collision of a soliton with the vertex there exists a time lapse during which the dynamics 
can be described as the scattering of three splitted solitary waves, one reflected on the same edge 
where the soliton was resident asymptotically in the past, and two transmitted solitary waves on 
the other edges. The amplitudes of the reflected and transmitted solitary waves are given by the 
scattering matrix of the linear dynamics on the graph. This behaviour takes place with a small error 
along a time scale of the order In v after the collision, where v is the asymptotic velocity of impinging 
soliton. On this time scale. Figure 4 is a realistic approximation of the process. 

t = t large 





Figure 4. Scattering of fast soliton on a Y-junction 
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The results in [22] and here discussed are inspired by the analogous analysis for NLS equation on the 
line with a repulsive 6 potential in the paper [52]. It should be said that a graph with two edges is 
equivalent to a line with a point interaction, so the treatment in [22] where several types of vertices 
are considered (Kirchhoff and repulsive 6 and 6', but the analysis could be extended to more general 
boundary conditions) shows how to generalize the results of [52] to other point interactions on the 
line and star graphs. To simplify the exposition, let us consider the following setting 

• Cubic NLS (this is essential: see later) 

• Kirchhoff vertex Ha = H (more general boundary conditions are allowed) 

• Initial state {v 3> 1; xq > v^~^, with < 6 < 1; x ^ smooth cutoff of the tail at the vertex) 

(2.21) ^o(a;) = iV2xix)e~'^'' cosh"i(x - Xq), 0, 0) 

We are interested in the evolution "^t of this initial condition. To this end, we will find an approximate 
solution of the equation: 

(2.22) ^t = e-'^'^Ho + i f e-'^(*-^)|^,p^, rfs . 

Jo 

The dynamics can be divided into three phases. 

The first phase is the approach to the vertex in the time interval t G [0, ti], where ti = Xq/v — v^^ . 
In this phase the incoming (quasi) soliton moves from xq to Xq — vti = v^~^ and ends the run at 
a distance of order v^~^ from the vertex. During this phase only a small tail of the pulse touches 
the vertex, the solution \I't behaves much as the solitary solution of the NLS in M and it remains 
supported on the edge ei with an exponentially small error. Choosing as the approximating function 
^t{x) = {(f)xo,-v{x,t),0,0) the following estimate (in L^ norm: control of masses) holds true 

Lemma 2.2. For any t G [0,ti] , ||\l>i - <l>t|| < Ce""'"* /or t; > 1 and 6 G (0,1). 

The proof consists in the accurate use of the well known and already cited Strichartz estimates (an 
ubiquitous tool in the study of nonlinear dispersive equations: see ^T]) to control distance between 
the unperturbed NLS flow and the NLS flow on the graph. 

The second phase is the interaction phase, when the "body" of the soliton crosses the vertex. It 
occurs during the time t G [ti, t2], where ^2 = xq/v + v^^]. The time interval ^2 — ^i is small (of order 
v~^, < 6 < 1), and the effect of the nonlinear term is demonstrably small. The soliton is fast and 
the pulse travels for a large distance (of order v(t2 — ti) = v^~^, < S < 1). This allows the linear 
dynamics being described by using a scattering approximation. Let ^ and M be the transmission 
and reflection coefficients of the Kirchhoff interaction H, ^ = ^ and M = —^^-^ ■ This means 
that the function '^{k,xi,X2,X3) = (e~*'^^^ +^6*^^^^ S/'e^^^'^, a^^ikx-i^ satisfies the Kirchhoff boundary 
conditions, and it is a solution of if \I' = A;^\l/ in distributional sense. The approximating function is 



chosen then as (recall 2.20) 



"^f = (0xo,-.W +=^0-x-o,.(i), '^^-xo-.W, '^^-.'o.^W) 

The choice of this reference approximate dynamics is explained in Figure 5. Reflected and trasmit- 
ted contributions are represented as tails (with ^ and S/' factors) of travelling solitons on "ghost" 
halflines. This trick is useful and has the further advantage to introduce fictitious lines where some 
representations and known properties of solitary waves are at disposal (see [22] for details). In any 
case, the true solution is compared with this approximate solution with an error small in an inverse 
power of velocity: 

Lemma 2.3. For any t G [ti,t2] \^t - $f || < Cv'^^'^ forv'^X and 6 G (0,1). 

The proof is technical and makes use of various well chosen representation for the linear time 
evolution in the interaction phase, an integral equation for the difference between the true solution 
and the approximate one, and an iteration of Strichartz estimates for evaluating errors. 
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. ^y^ 




Figure 5. On the left side the state $'^ at time ti. On the right side the state $"^ at 
time ^2- Each edge of the graph is extended to a hne by ideally adding the half line 
(— oo, 0] represented by the dashed lines. 

Finally there is the post interaction phase, in the time interval t G [^2)^3]? with t^ = t2 + T\Y\.v where 
the free NLS dynamics dominates again; however, now the initial data are not asymptotic solitary 
waves, but waves with approximately soliton-like profiles and wrong amplitudes, due to the presence 
of scattering coefficients. Precisely, for t >t2 and x G M we define the functions 0*^ and cj)^'^^ by 
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The approximate solution is defined as 
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and the it satisfies the following lemma, which is the main result. 

Lemma 2.4. Fix T > 0, then for any time t G [^2,^2 + Tint"], there exists < r] < 1/2 such that 

ii^i-$r*ii <cv-'^ 

To prove the lemma one has to use in an essential way the fact that the cubic NLS equation is 
integrable on the line, and thanks to this it is possible to get large time behaviour of initial "lowered 
solitons" along the nonlinear evolution (see [52] which is the source of the idea). Again with due 
estimates of the errors these results can be translated on the graph, ending with the quoted result. 
We end with some remarks and comments. 



The linear Hamiltonians in 2.22 to which the theorem refers, when different of the Kirchhoff one, 
have to be rescaled in order to give a nontrivial scattering matrix in the regime of high velocity; for 
example in the case of a delta potential one has to set a — )■ f a . 

The estimates are in L^-norm. So the present analysis of soliton scattering is rigorous for what 
concerns mass transmission and reflection. On the contrary, the result proposed in [52j is in L°°- 
norm and the control is on the proflle of the outcoming pulses. 



In the time interval t2 < t < t2+Tlnv , the last lemma implies that 






\^\ + 0{v-'^) (j = 2,3) 



holds for a certain a > and where ^ is the transmission coefficient. So, in the limit of fast solitons, 
i.e. V — )■ 00, the ratio which defines the nonlinear scattering coefficient converges to the corresponding 
linear scattering coefficient. Analogously for refiection coefficient. 

The results discussed in this section are given and proved for a repulsive interaction (no negative 
eigenvalues) at the vertex. In the presence of eigenvalues a more refined analysis has been performed 
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in [53] for the NLS equation on the hne with an attractive S potential, giving resuhs quahtatively 
similar to the repulsive case. Notice that this is the case in which stable standing waves exist. 
In view of the previous remarks, the directions in which the scattering of solitons on star graphs 
could be extended with a certain amount of technical work, but without introducing new ideas are: 
a)the (straightforward) case of NLS on star graphs with more than three edges; the case in which 
the underlying linear quantum graph with 6 or 6' vertex has bound states; with some caution, the 
case of general selfadjoint boundary conditions. Nothing it is presently known about the possible 
existence of multisolitons, or about the collision of several solitons at the vertex. 
On the contrary, scattering of solitons on graphs having a less trivial topology is an open problem. 
In particular it could be interesting the study of propagation of solitons on trees, about which 
something has been said before in relation to BBM equation. Notice that among the main technical 
ingredients are Strichartz estimates, and every generalization needs their validity. Recent advances 
in this direction are given in [261 I2Z!- To give an idea of the interplay between nonlinearity and 
scattering on complex networks, let us briefly discuss the interesting paper [3] . The authors consider 
a (possible complex) graph with bounded edges where NLS dynamics is posed, and to this localized 
nonlinear network two external legs are attached where linear propagation occurs. After some general 
remarks, numerical results are given and discussed concerning the scattering of stationary wave 
OinC**^^ incoming from one of the external edge, entering the nonlinear network and outgoing from 
the second external edge: numerical experiments are done for a large spectral range of fc's and 
intensities Im = |amp . It turns out that also for relatively simple examples of networks, scattering 
is dominated by sets of sharp resonances of the underlying linear model; these tend to sensibly 
amplificate the effect of nonlinearity and they prevent to consider it as a small perturbation. This 
moreover gives rise to typical effects in nonlinear dynamics, such as multistability and hysteresis. 
The resonances correspond to longliving states captured in the networks, and they are a result of 
the non trivial topology of the graph. In this sense the effect should be considered as a general 
phenomenon which could be expected in real experimental settings involving nonlinear networks. 
Another important problem concerns the possibility of extending the timescale of validity of approx- 
imation by the solitary outgoing waves. In a different model (scattering of two solitons on the line) 
in [5l], some considerations are given about the possibility of longer timescales of outgoing soliton 
approximation depending on the initial data and external potential, but it is unclear whether similar 
considerations can be applied to the present case. 

As regards the limitation to cubic NLS, the fundamental asymptotics proved in [32] and used to get 
the analogous result on graphs depend on the integrability of cubic NLS. One can conjecture that for 
nonlinearities close to integrable the outgoing waves are close to solitons over timescale similar to the 
above ones. In this respect it appear as interesting the recent result of Perelman on the asymptotics 
of colliding solitons on the line ( [55] ) . 

As a final remark, let us notice that the two main dynamical features here described and about 
which a description, if partial, has been achieved, i.e. standing waves and their neighborhood on one 
hand, and scattering of fast solitons on the other hand, correspond to states and regimes far apart in 
the energy space of the system. Many interesting physical phenomena experimentally well tested are 
probably in a different or intermediate region. In particular capture and generation of solitons are 
not yet understood from a rigorous point of view, and they represent a challenge to mathematical 
methods and theoretical interpretation where perhaps the simple but non trivial model of NLS 
equation on graphs could offer some insight. 
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